HD 137510: AN OASIS IN THE BROWN DWARF DESERT 
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ABSTRACT 

Since the beginning of precise Doppler surveys, which have had stunning success in detecting 
extrasolar planetary companions, one surprising enigma has emerged: the relative paucity of 
spectroscopic binaries where the secondary mass lies in between the stellar and planetary mass 
' regime. This gap in the mass function for close-in (a < 3 — 4 AU) companions to solar-type 

stars is generally referred to as the "Brown Dwarf Desert" . Here we report the detection of a 
companion to HD 137510 (GOIV), with a minimum mass of 26 Mj up it r, moving in an eccentric 
orbit (e = 0.4) with a period of 798 days and an orbital semimajor axis of 1.85 AU. The detection 
is based on precise differential radial velocity (RV) data obtained by the McDonald Observatory 
and Thiiringer Landessternwarte Tautenburg planet search programs. 

Subject headings: stars: late-type — stars: low mass — planetary systems — techniques: radial velocities 
1. Introduction 



1 Based on observations made at McDonald Observatory 
and the Thiiringer Landessternwarte Tautenburg It was the common expectation that, with the 

2 current address: Department of Astronomy, University increase in the precision of RV measurements, the 

of Michigan, 830 Dennison Ann Arbor MI 48109, USA detected secondary masses of spectroscopic bina- 

3 current address: Astrophysikalisches Institut, Univer- rieg will gra d u ally decrease, first into the brown 

sitat Jena, Schillcrgasschen 2-3, D-07745 Jena, Germany , r ,-. , „„ , , % , 

dwarf mass range (M < 80 Mj up ) and eventu- 
ally into the planetary mass regime of roughly 
M < 13Mj up (below the deuterium burning limit). 
Quite the opposite turned out to happen. The pi- 
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oneering long-term Dopplcr survey of Campbell, 
Walker, & Yang (1988), and later also in the 
southern hemisphere by Murdoch, Hearnshaw, & 
Clark (1993), did not reveal any brown dwarf com- 
panions to nearby solar-type stars, despite their 
sufficient sensitivity. There appeared to be an 
abrupt drop in the frequency of companions at the 
hydrogen burning limit of w 0.075 M Q . The first 
very low-mass companion (msini = 11 Mj up ) de- 
tected by the RV technique is the companion to 
HD 114762 (Latham et al. 1989). The discoverers 
classified the object as a brown dwarf compan- 
ion, but according to Cochran, Hatzes, & Han- 
cock (1991) this companion might be in fact a 
late M dwarf due to a possible low value of the 
inclination angle i. In the years after the discov- 
ery of the first unambiguous extrasolar planet or- 
biting a main-sequence star by Mayor & Queloz 
(1995), a total of more than a 100 extrasolar 
planets have been accumulated by several precise 
Doppler searches all around the globe. Still, with 
a few exceptions (HD 127506: Mayor, Queloz, 
& Udry 1998; HD 29587 & HD 140913: Mazeh, 
Latham, & Stefanik 1996; HD 168443 (second 
companion): Marcy & Butler 2001; HD 184860: 
Vogt ct al. 2002), the part of the mass function 
for spectroscopic secondaries, in between the min- 
imum hydrogen burning mass and the deuterium 
burning limit of M « 13 Mj up , remained surpris- 
ingly sparsely populated. A notable exception is 
HD 10697, where Vogt et al. (2000) detected a 
companion with msini = 6.35 Mj up , while Zuckcr 
& Mazeh (2000) determined a true mass for the 
secondary of 40 Mj up after combining the spectro- 
scopic solution with Hipparcos astrometry. Udry 
et al. (2002) found 4 companions with minimum 
masses close to the planet/brown dwarf border, 
which could be either "superplanets" or real brown 
dwarfs. However, most of the candidate brown 
dwarf secondaries considered in the Halbwachs et 
al. (2000) study turned out to be stellar mass com- 
panions, again after combination with Hipparcos 
astrometry. 

This suggests the existence of two distinctive bi- 
nary formation processes for certain mass-ratios. 
The gap in the mass function between the stel- 
lar and planetary mass range for close binaries is 
generally dubbed the "Brown Dwarf Desert" (this 
desert does not appear to exist at wide separa- 
tions, e.g. Gizis et al. 2001: Ncuhauser & Guen- 



ther 2004). The gap is particularily obvious as 
a sudden drop at the deuterium burning limit in 
the substellar mass function for close companions. 
While planets are relatively common as compan- 
ions to solar-type stars, brown dwarfs are rare. 
This fact is especially surprising since brown dwarf 
companions are much easier to detect, even for 
lower precision surveys, due to the larger RV am- 
plitudes. Moreover, these objects cover a more 
than 4 times larger mass range than planetary 
companions, still their detected numbers per mass 
bin is extremely low compared to lower and higher 
mass companions. 

In this work we present the discovery of another 
possible "oasis" in this desert, a brown dwarf or- 
biting the G-type star HD 137510. 

2. Stellar characteristics of HD 137510 

HD 137510 (= HR 5740 = HIP 75535) is a 
bright (V = 6.3) G-type star at a distance of 
41.75 ± 1.7 pc, according to the Hipparcos par- 
allax of 23.95 ± 0.94 mas (ESA 1997). In order 
to derive the stellar parameters for the primary 
we used the (B-V)-T e ff calibration of Alonso, Ar- 
ribas, & Martinez-Roger (1996), compared My 
and T c ff with stellar evolution calculations (e.g. 
Rcddy et al. 2003), and performed a model atmo- 
sphere analysis of a high resolution, high signal- 
to-noise spectrum of HD 137510. We find that the 
star has a mass of M = 1.3 ± 0.1 M Q , an effec- 
tive temperature T c ff = 5896 ± 57 K, a log g of 
4.0 ± 0.1, a radius R of 1.9 ± 0.2 R Q and an age of 
3.4 ±0.8 Gyrs. An iron abundance analysis yields 
[Fe/H] = 0.16±0.07 dex. The object is thus metal 
rich compared to the Sun. We also find substan- 
tial rotational broadening of the spectral lines, the 
value of wsini can be constrained in the range 4.5 
to 7.5 krns" 1 . These parameters are consistent 
with (or very close to) those determined by Lebrc 
ct al. (1999) and do Nascimento et al. (2000). It 
thus appears that HD 137510 was a metal rich 
late F-type star which has now started its evolu- 
tion away from the main sequence. Interestingly, 
Smith, Cunha, & Lazzaro (2001) have noted that 
HD 137510 possibly underwent chemically frac- 
tioned accretion and represents a good candidate 
for having a close-in giant planetary companion. 
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3. Radial velocity data and Keplerian or- 
bital solution 

We carried out observations of HD 137510 as 
part of the ongoing Doppler surveys for extraso- 
lar planets using the Harlan J. Smith 2.7 m tele- 
scope at McDonald Observatory (McD) in West 
Texas and the 2.0 m Alfred Jcnsch Telescope at 
Thiiringer Landessternwarte Tautenburg (TLS) in 
Germany. In September 1998 a first spectrum of 
HD 137510 was taken at McD, but RV monitoring 
was discontinued for almost three years. In May 
2001 observations of HD 137510 were started at 
TLS, while at the same time monitoring was re- 
sumed at McD. For both programs we use a tem- 
perature stabilized molecular iodine vapor (I 2 ) cell 
as the velocity metric. To obtain the high preci- 
sion differential RV measurements we employ spe- 
cialized data modeling techniques, which include 
the detailed reconstruction of the instrumental 
profile of the respective spectrometers at the time 
of observation (Valenti, Butler, & Marcy 1995; 
Butler et al. 1996; Endl, Kiirster, & Els 2000). 

Fig. 1 displays our RV data as a function of time 
and with the best-fit Keplerian orbital solution su- 
perimposed. We employ the program Gaussfit 
(McArthur, Jefferys, McCartney 1994), which uses 
a robust estimation method to find orbital so- 
lutions for each data set (McD and TLS) inde- 
pendently, as well as a combined solution, where 
the arbitrary velocity zero points (7 velocities) of 
both data sets are included as fit parameters. Ta- 
ble 1 summarizes the parameters we determined 
for the orbit. The resulting mass function of this 
binary yields a minimum mass for the companion 
of 26 ± 1.4 Mj up . Table 2 lists our Doppler mea- 
surements for HD 137510. 

The residual scatter around this orbital fit is 
unusually high, McD: a = 13.3 ms" 1 and TLS: 
(T = 21.4 ms -1 , roughly a factor of 2 worse than 
the routine RV precision for each telescope system. 
It also exceeds the average internal measurement 
errors of both data sets (McD: 8.9 ms" 1 , TLS: 
12.0 ms -1 ). Stellar activity can introduce addi- 
tional noise in RV data at that level of measure- 
ment precision (e.g. Saar & Donahue 1997). For 
HD 137510 we can rule out this possibility, based 
on the general low Ca II H&K emission of this 
star. The McD spectra, for which we determined 
the RVs, also include the Ca II H&K lines, two well 



known activity indicators. We find a McD S-indcx 
SmcD = 0.157±0.011 (see section 2.3 in Paulson et 
al. 2002 for a detailed description of how we obtain 
the S-index for the Ca II K line core emission). 
This is a comparable emission level as the well 
known inactive star r Ceti (SmcD = 0.166 ±0.008) 
and shows that HD 137510 is not an active star. 
Furthermore, the star does not show up as an X- 
ray source in the ROSAT all-sky survey results 
(Hiinsch, Schmitt, & Voges 1998). One possible 
explanation for the higher scatter is the substan- 
tial rotational broadening of the spectral lines in 
the case of HD 137510 (vsini = 4.5-7.5 kms -1 ). 
The larger scatter for the TLS data likely reflects 
the differences both in telescope aperture (TLS: 
2.0 m ; McD: 2.7 m) as well as in the relative see- 
ing conditions at these two sites. Another possible 
explanation is the presence of a second low-level 
Keplerian signal. We searched in the RV residuals 
for periodic signals using the Lomb-Scargle peri- 
odogram (Lomb 1976, Scargle 1982). While the 
two separate periodograms of the individual sets 
of residuals do not reveal a significant signal, we 
do find a strong signal with a period of 227 days 
in the combined residuals. This warrants further 
intensive monitoring of HD 137510 to confirm a 
possible planetary companion at a separation of 
« 0.8 AU (which would make this system very 
interesting from a dynamical standpoint). 

4. Discussion 

Due to the general sparseness of brown dwarf 
RV companions we have to look for additional 
clues in order to rule out a stellar mass companion 
producing the detected Keplerian signal: 

Simply based on the probabilty of a random 
distribution of the inclination angle i, we already 
have a 90% probability that i is larger than 25°, 
which translates into a true mass of < 61 Mj up , 
still below the hydrogen burning limit. 

We measure a vsini of 4.5 to 7.5 kms -1 for 
HD 137510, which is a typical value for a G0IV 
star. The mean vsini value of the F8 - G2 sub- 
giants in the de Medeiros & Mayor (1999) cata- 
logue is 6.5 kms -1 (excluding one fast rotating 
F8IV star the mean value drops to 5.7 kms -1 ). 
This shows that the value for HD 137510 is not 
atypically low, which means that we probably do 
not have a near pole-on view on the primary star. 
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Assuming that the orbital plane is roughly perpen- 
dicular to the rotation axis of the star we, again, 
do not seem to face a low i scenario in this case. 

McAlister et al. (1987) obtained a negative re- 
sult for this star in their speckle survey for duplic- 
ity among Bright Star Catalogue stars using the 
CFHT. Granted, a positive detection would have 
been a difficult one, close to the diffraction limit 
of the telescope (« 0.04 arcsec). 

HD 137510 has no flag in the multiple/binary 
annex of the Hipparcos catalogue, showing that 
no suspicious astrometric variation was discov- 
ered. According to Lindegren et al. (1997) a 
stochastic solution is chosen for objects where no 
acceptable single or double star solution can be 
found and additional noise, typically ranging from 
3 to 30 mas, needs to be added until an acceptable 
fit is found. If we regard the 3 mas as the lower 
limit of an astrometric pertubation (angular size 
of the semimajor axis) which would have shown up 
in the analysis of the Hipparcos data of a bright 
target like HD 137510, we can exclude inclination 
angles of i < 16° and a true mass of the compan- 
ion of > 94 Mjup, using the constraint of Pourbaix 
& Jorissen (2000). 

5. Conclusion 

We have presented observational evidence 
for the substcllar nature of the companion to 
HD 137510. The question, why this star consti- 
tutes one of the very rare oases in the brown dwarf 
desert, remains unanswered. However, due to the 
well established parameters of the primary and the 
orbit, this new system might help to understand 
binary formation processes for these mass-ratios 
at small orbital separations. 

We are grateful to the McDonald Observatory 
Timc-Allocation-Committee for generous alloca- 
tion of observing time. This material is based 
upon work supported by the National Aeronautics 
and Space Administration under Grant NAG5- 
9227 issued through the Office of Space Science, 
and by National Science Foundation Grant AST- 
9808980. We thank our referee Adam Burgasser 
for his helpful comments. 
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Table 1 

Orbital parameters for the companion to HD 137510. Solutions are given for both data 

SETS SEPARATELY, AS WELL AS THE COMBINED SOLUTION. 





McD 


TLS 


McD+TLS 


Period [days] 


799.0 ± 1.4 


804.5 ± 23.2 


798.2 ± 1.4 


Tperiastron [JD] 


2452583.2 ± 3.2 


2452572.3 ± 9.5 


2452582.01 ± 2.6 


Kims" 1 ] 


531.7 ± 6.9 


508.9 ± 22.5 


531.6 ± 5.3 


e 


0.40 ± 0.008 


0.397 ± 0.03 


0.402 ± 0.008 


UJ 


30.9 ± 1.4 


28.7 ± 3.1 


30.8 ± 1.2 


msini [Mj up ] 


26 ± 1.4 


25 ± 1.4 


26 ± 1.4 


a [AU] 


1.85 ± 0.05 


1.86 ± 0.05 


1.85 ± 0.05 


rms [m s _1 ] 


13.1 


21.0 


McD:13.3 TLS:21.4 
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Table 2 

Differential radial velocities for HD 137510 (note that both data sets have arbitrary 

AND THUS DIFFERENT VELOCITY ZERO-POINTS) . 



JD [2400000+] dRV [ms -1 ] error [ms *] Survey 



51067.6406 


165.7 


7.5 


McD 


52037.7836 


-209.7 


7.5 


McD 


52040.8490 


-182.1 


7.1 


McD 


52114.7432 


-154.0 


7.5 


McD 


52141.6609 


-133.2 


8.2 


McD 


52142.6541 


-151.1 


8.3 


McD 


52144.7022 


-133.4 


17.4 


McD 


52145.6828 


-123.6 


8.3 


McD 


52331.9521 


101.9 


7.7 


McD 


c n q ir Q no 1 n 


134.5 


15.2 


IV r „ y\ 
IviCU 


52357.9625 


151.6 


6.6 


McD 


52388.9414 


265.5 


9.1 


McD 


52451.7614 


473.5 


7.4 


McD 


52472.7426 


577.6 


7.6 


McD 


52492.6720 


665.2 


7.4 


McD 


52576.5355 


822.8 


9.3 


McD 


52689.0031 


62.4 


9.4 


McD 


52689.0124 


53.4 


8.9 


McD 


52741.8983 


-124.1 


11.9 


McD 


52741.9064 


-131.8 


9.4 


McD 


52743.9543 


-119.9 


7.7 


McD 


52804.8642 


-184.4 


7.3 


McD 


52804.8749 


-183.0 


7.5 


McD 


52807.6845 


-191.4 


7.7 


McD 


52807.6914 


-190.1 


7.5 


McD 


52839.7200 


-186.1 


8.3 


McD 


52895.5874 


-174.6 


8.5 


McD 


52895.6010 


-168.5 


8.0 


McD 


52896.6072 


-179.5 


7.4 


McD 


52931.5547 


-176.0 


8.2 


McD 


52931.5686 


-180.9 


8.4 


McD 


53018.9949 


-83.5 


9.2 


McD 


53018.0065 


-85.7 


8.8 


McD 


53038.9840 


-27.5 


8.7 


McD 


52039.5781 


-220.5 


7.1 


TLS 


52040.5039 


-233.4 


6.4 


TLS 


52042.4922 


-224.5 


6.4 


TLS 


52066.4219 


-239.6 


7.6 


TLS 


52090.4375 


-186.8 


6.8 


TLS 


52213.2188 


-81.1 


18.7 


TLS 


52280.6875 


-31.6 


11.0 


TLS 


52334.5703 


77.4 


9.5 


TLS 


52365.5430 


162.0 


6.3 


TLS 


52366.5469 


173.5 


8.6 


TLS 


52394.5391 


280.6 


9.7 


TLS 


52448.4219 


445.9 


5.7 


TLS 


52452.3984 


471.8 


7.7 


TLS 


52478.4805 


529.5 


19.8 


TLS 


52483.4062 


573.5 


16.1 


TLS 


52506.3320 


656.7 


15.7 


TLS 


52516.3125 


754.5 


14.3 


TLS 


52517.2969 


722.5 


28.6 


TLS 


52519.3086 


762.8 


12.5 


TLS 


52684.6445 


15.5 


14.0 


TLS 


52745.5664 


-156.9 


11.0 


TLS 


52754.5312 


-174.1 


8.5 


TLS 


52804.4180 


-217.1 


10.3 


TLS 


52831.3984 


-187.3 


7.6 


TLS 
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Fig. 1.— RV data for HD 137510 (GOIV), from 
McDonald Observatory (McD) and Tautenburg 
Observatory (TLS). The best-fit combined Kep- 
lerian orbital solution (see table 1 for the orbital 
parameters) is plotted as solid line. The lower 
panel shows the residuals after subtraction of this 
orbit. 
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Table 2 — Continued 



JD [2400000+] dRV [ms -1 ] error [ms -1 ] Survey 



52833.3984 


-224.4 


8.5 


TLS 


52834.4102 


-193.4 


10.5 


TLS 


52835.3984 


-192.3 


8.4 


TLS 


52836.4766 


-195.8 


10.2 


TLS 


52838.4766 


-198.3 


14.6 


TLS 


52839.4141 


-187.7 


9.9 


TLS 


52840.4023 


-190.5 


10.8 


TLS 


52841.3516 


-225.3 


12.3 


TLS 


52861.3750 


-185.1 


13.2 


TLS 


52863.3516 


-237.1 


17.1 


TLS 


52872.3398 


-211.5 


18.7 


TLS 


52877.3672 


-193.1 


10.3 


TLS 


52878.3594 


-189.0 


13.0 


TLS 


52930.2266 


-195.6 


9.8 


TLS 


52931.2227 


-222.3 


7.2 


TLS 


52952.1992 


-135.6 


11.2 


TLS 


52956.1992 


-164.3 


16.1 


TLS 


52982.7031 


-144.2 


12.5 


TLS 



9 



